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ABSTRACT 

Context. Tracing the evolution of debris discs is essential for our understanding of planetary system architectures. While the evolution 
of their inner discs has been recently studied with the Spitzer Space Telescope at mid- to far-infrared wavelengths, the outer discs are 
best characterised by sensitive millimetre observations. 

Aims. The goal of our study is to understand the evolution timescale of circumstellar debris discs, and the physical mechanisms 
responsible for such evolution around solar-type stars. In addition, we perform a detailed characterisation of the detected debris discs. 
Methods. Two deep surveys of circumstellar discs around solar-type stars at different ages were carried out at 350 /jm with the CSO 
and at 1.2 mm with the IRAM 30-m telescope. The dust disc masses were computed from the millimetre emission, where the discs 
are optically thin. Theoretically, the mass of the disc is expected to decrease with time. In order to test this hypothesis, we performed 
the generalised Kendall's tau correlation and three different two-sample tests. A characterisation of the detected debris discs has been 
obtained by computing the collision and Poynting-Robertson timescales and by modelling the spectral energy distribution. 
Results. The Kendall's tau correlation yields a probability of 76% that the mass of debris discs and their age are correlated. Similarly, 
the three two-sample tests give a probability between 70 and 83% that younger and older debris systems belong to different parent 
populations in terms of dust mass. We detected submillimetre/millimetre emission from six debris discs, enabling a detailed SED 
modelling. 

Conclusions. Our results on the correlation and evolution of dust mass as a function of age are conditioned by the sensitivity limit of 
our survey. Deeper millimetre observations are needed to confirm the evolution of debris material around solar-like stars. In the case 
of the detected discs, the comparison between collision and Poynting-Robertson timescales supports the hypothesis that these discs 
are collision dominated. All detected debris disc systems show the inner part evacuated from small micron-sized grains. 

Key words, circumstellar matter - planetary systems: formation - stars: late-type - Kuiper Belt 



T""! ■ 1. Introduction discs, the dust grains are removed from the optically thin disc by 

rN radiation pressure, Poynting-Robertson drag or dust sublimation 

H ; As part of the star formation process in molecular cloud cores, timescales shorter t han th e stellar pre-main sequence lifetime 

^ circumstellar discs form conserving the initial angular mo- (Backman & Paresce 1993", Meyer et al. lOOT). Nevertheless, 

mentum. These so-called primordial discs have typical masses ^n infrared excess above the stellar photosphere has been dis- 

(10 -10- Mo) and sizes comparable to that expected for the covered around more than 300 pre- and main-sequence stars, 

primitive solar nebula and they provide the environment and jhese stars are th ought to host a remnant debris disc. Backman 

mater ial fro m which pl anets a re expected to form (Beckwith & Paresce ( fT993T l suggested that the infrared excess of these ob- 

et al. |2000| Natta et al. |2007|»- The first stages of planet for- je^ts comes from small dust particles which ai-e the products of 

mation occur close to the disc mid-plane where sub-micron collisions of larger bodies. Thus, a debris disc should consist of 

grains grow to millimetre sizes and settle (e.g. Beckwith et ^ 'second generation' of small dust paiticles, larger bodies and, 

al. ,2000 ). Larger planetesimals may successively form by colli- possibly planets 

sions between grains and through gravitational instabilities (e.g. xhe formation of large bodies is thought to be f aster in the 

IJohansen et al.2Q07j Henningi2008il. The primordial gas in cir- inner part of the disc (jKenyon & Bromley 2004a). As a 



con- 



cu mstell ar discs is th ought t o last for 1-10 Myr (e.g. Haisch et sequence, new small dust particles will be produced by colli 
al- 12001 1 Lawson et al.|2004|l. After the gas is dispersed from the sions first in the inner region of the disc and, later, further out 



(Kenyon & Bromley 2004a I. During this phase, small dust par- 



Send o^nw? re^wei-f^ to: V. Roccatagliata: roccata@mpia.de tides are continuously produced in the disc. In this perspec- 

* This work is based on observations made with the IRAM (Institut tive, transitional objects are of particular interest. These may 

de Raclioastonomie MilUmetrique) 30-m telescope and the CSO represent a transitional phase between the primordial optically 

(Caltech Submillimetre Observatory) 10-m telescope. IRAM is sup- thick disc phase and the debris optically thin disc phase. Their 
ported by INSU/CNRS (France), MPG (Germany), and IGN (Spain). 
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infrared excess suggests the presence of an inner hole in the 
disc where the small grain population, which emits in the near- 
IR, is missing. Some of them are still accreting material onto 
the central star, while others show intense collisional activity in 
their inner part, possibly induced by the formation of planets 
jCurrie et al. 20071 [Keiiy on & Bromley 2004b| l. 
While the dissipation of debris discs takes place by radiation 
pressure and/or Poynting-Robertson drag and the total dust mass 
decreases with time, at early stages the stellar wind drag can 
be an important dust removal mechanism. Depending on which 
of these processes is the main driver of dust removal, we ex- 
pect a power-law dependence with a different index (Dominik & 
Decin (2003 ), Wyatt et al. ( 2007a ,b), Lohne et al. (20 08 )) and 
a different dust distribution if a planetesimal belt at a given dis- 
tance is creating dust of just one size (Wyatt 2005i Spangler 
et al. ( 1200 11 1 measured the fractional luminosity (LiR/Lstar) of 
pre- and main sequence stars with ISO. They found that the 
fractional luminosity, which is proportional to the dust present 
in the disc, decreases with time as f While the analysis of 
Spangler et al. (2001) included also pre-main sequence stars, 
Habing et al. (2001 ) analysed the ISO and IRAS data of nearby 
main-sequence stars alone: they concluded that most of the stars 
reach the main-sequence surrounded by a disc which decays 
in about 400 Myr Decin et al. (I2003I I reviewed the different 
studies on the time dependency of Vega-like systems done with 
ISO/IRAS. They concluded that the observations only indicated 
a large spread in fractional luminosity for stars at most ages, 
that the maximum excess was independent of time and that the 
debris disc phenomenon was more common in younger stars. 
An analysis of the near to far infrared excess, carried out with 
Spitzer, led to the finding of a general decrease of the fractional 
luminosity with time as r^, starting later at longer wavelengths 
jSu et al. 2006i l. This suggested a disc clearing more efficiently 
in the inner parts dSu et al. 20061 ). Deviations from these trends 
have been associated with a recent or still ongoing collisional 
cascade that produces a small dust grain population which is 
rapidly removed by radiation pressure or by the action of the 
Poynting-Robertson drag (S ong et al. 2005 Wyatt et al. l2007bl) . 
A complementary picture of the evolution of the outer part of the 
debris discs is provided by (sub)millimetre observations. These 
show a decline in the dust mass as a function of age and signif- 
icant evolution of the circumstellar dust mass within the first 10 
Myr and between 10 Myr and a few Gyr (e.g. ILiu et al. 20041 
Carpenter et al. 12005] (hereafter C05)). 

In this paper we present the results of two deep sub-millimetre 
and millimetre surveys of circumstellar discs around solar-type 
stars with ages between 3 Myr and 3 Gyr. The characteristics 
of our sample are presented in §2.1. The observations carried 
out at 350 jum and at 1.2 mm are reported in section 2.2 and 
section 2.3. The results from our two surveys are presented in 
section 2.4. The discussion of the dust disc mass evolution as a 
function of the age of the systems, compared with the history of 
our solar system, is presented in section 3. A detailed analysis of 
the debris disc properties by modelling their spectral energy dis- 
tribution (SED) is provided in section 4 and finally the summary 
is presented in section 5. 



2. Observations 

2.1. Sample 

Nearby sources (d<150 pc) were selected from The Formation 
and Evolution of Planetary Systems (FEPS) 5/?/fzer Legacy pro- 
gram (Meyer et al. l2006l l. The FEPS sample contains 314 stars 



with stellar masses between ~0.5 and 2 Mq, ages spanning the 
range from 3 Myr to 3 Gyr and spectral types between K7 and 
F5. In addition, 14 objects were initially included in the FEPS 
sample because of an infrared excess previously detected with 
IRAS/ISO. Since these sources can introduce a bias in the FEPS 
sample with respect to the presence of a disc, they are not taken 
into account in any statistical analysis in this paper 
The characteristics of the entire star sample have been exten- 
sively analysed by the FEPS team (e.g. Meyer et al. 120061 
Mamajek & Hillenbrand |2008j). The stellar ages are based on 
pre-main sequence tracks for stars younger than 10 Myr, X-ray 
activity and the strength of the Call H and K emission for nearby 
solar-type stars (<50 pc), and the association of stars with clus- 
ters or star-forming regions of known age. A new calibration 
presented in Mamajek & Hillenbrand ( 120081) has been applied. 
Final ages will be presented by Hillenbrand et al. (in prepara- 
tion, see also Mamajek & Hillenbrand 2008). The distances of 
the sources were determined by the FEPS team on the basis of 
Hipparcos parallaxes for nearby stars and kinematic distances 
for stars associated with young moving groups and associations. 
Our total sample contains 141 sources of which: 

- 16 targets were observed at 350 /urn with the CSO {Caltech 
Submillimetre Observatory) 10-m telescope (see col. 6 Table[TJ 

- 40 targets were observed at 1.2 mm with the IRAM (Institut 
de Radioastronomie Millimetrique) 30-m telescope (see col. 7 
Table □ 

- 121 targets are taken from the SEST (Swedish-ESO 
Submillimetre Telescope) and OVRO (Owens Valley Radio 
Observatory) surveys performed by Carpenter at. al (120051 ). of 
which 17 have been re-observed during our surveys. 

The sample includes 15 sources younger than 10 Myr, 16 
intermediate-age systems (10-20 Myr) and 110 evolved systems 
(> 20 Myr). In section 3 we only consider the evolved systems in 
order to search for a correlation between the mass of the debris 
disc and the age of the system. In particular, 104 out of the 110 
evolved systems are analysed excluding the sources not compati- 
ble with the FEPS sample selection (see above). The 104 sources 
analysed here include 45% of the FEPS sample older than 20 
Myr and in particular: 35% of the systems with no infrared ex- 
cess (until 70 yum) and 89% with excess. 



2.2. CSO observations at 350 yum 

Observations of the sources listed in Table 1 were carried out at 
the CSO between 17 and 21 April 2005, using the Submillimetre 
High Angular Resolution Camera II (SHARC-II) with the 350 
fim filter The full width at half maximum (FWHM) for the 
10.4 m telescope beam profile is 8.5" at this wavelength. The 
median zenith opacity at 225 GHz was 0.04. The typical on- 
source integration times varied from 40 minutes to 2 hours 
depending on source elevation and zenith opacity, although a 
deeper 3.5 hour integration was obtained for HD 107146 (see 
Corder et al. ,2009) . Pointing and instrumental flux calibra- 
tors were selected from the CSO SHARC-II list of calibra- 
tors. Data imaging and reduction were carried out using the 
Comprehensive Reduction Utility for SHARC-II. The errors as- 
sociated with each detection are the quadratic sum of the instru- 
mental error (between 5.3 and 45 mJy) and the calibration errors 
(between 12% and 22% of the flux). Instrumental errors are com- 
puted based on the rms computed in a sky annulus, and propa- 
gated over the photometric aperture. Calibration errors were de- 
termined based on repeatability of the flux calibrators over the 
entire observing run. 



Roccatagliata et al.: Long- wavelength observations of debris discs around sun-like stars 



3 



2.3. IRAM observations at 1.2 mm 

Millimetre continuum observations of 40 sources at 1.2 mm 
were carried out between May and December 2005 (Projects 
num. 013-05, 141-05) and between January and April 2008 
(Project num. 185-07) at the IRAM 30-m telescope at Pico 
Veleta with MAMB02, the 117 channel IRAM bolometer The 
FWHM of the 30-m telescope beam is 11" at 1.2 mm. The me- 
dian zenith opacity at 230 GHz was 0.5 during the first run be- 
tween May and November 2005, and 0.3 during the second run 
carried out in December 2005 and during the last run in 2008. 
Pointing and calibrator sources were selected from a standard list 
provided by IRAM. The observations were carried out in on-off 
integration mode using a standard calibration. An rms of 0.7 mJy 
was obtained in one hour (blocks of three times 20 minutes) of 
on-off integration. The weather conditions were stable during 
most of the nights and the calibrators were observed within a 
few hours from the science targets. The data were reduced with 
the data reduction package MOPSIC (R. Zylka 2006). 
Observations of the same target at different dates were coadded 
to increase the signal to noise. The uncertainties associated with 
each flux are the quadratic sum of the statistical error and the 
calibration error. The statistical error is the rms of the flux over 
the entire telescope field of view, while the calibrator error is 
based on the agreement between the predicted calibrator fluxes 
and those measured over 48 hours before and after the obser- 
vations. We obtained statistical errors between 0.8 and 2.8 mJy 
and calibration errors between 11% and 16% in flux. The final 
uncertainties are dominated by the statistical error in the IRAM 
observations. 



2.4. Detections and upper limits 

The fluxes (in the case of detections) and the 3cr upper limits 
(in the case of non-detections) measured at 350 yum and 1.2 mm 
are listed in Table [1] The detection errors and the upper limits 
include only the statistical error. The calibration uncertainty is 
reported in the note to Table [T] The 3cr upper limits were com- 
puted as 3 times the statistical error described above. 31 sources 
have been observed with IRAM, 7 sources only with CSO and 
9 sources with both telescopes. We detected 1 1 sources in total: 
5 sources with IRAM, 8 with the CSO. Two of the sources were 
detected with both IRAM and CSO. Five of them are primordial 
discs younger than 10 Myr, and six are debris discs. All sources 
were unresolved with the exception of HD 107146 which had a 
deconvolved size at 350 fim of 8'.'9 x 8'.'2 with an uncertainty of 
about ^.'6 per axis (Corder et al. i2009l ). 



3. Disc dust mass 

The main goal of this project is to analyse the temporal evolution 
of circumstellar debris discs and the mechanism of their dissipa- 
tion around solar- type stars. 

The dust mass in debris discs is expected to decrease with time. 
Depending on which process (collisions or Poynting-Robertson 
drag) is the main driver of dust removal, we expect a different 
power-law dependence of the dust mass as a function of the age 
of the system (Dominik & Decin 12003]) . 

In this analysis we include the dust masses computed from the 
millimetre observations of nearby solar-type stars obtained with 
IRAM presented here and with SEST and OVRO (from C05). 
The disc is assumed to be isothermal and optically thin (r << 1) 



at millimetre wavelengths. The dust disc mass can then be de- 
rived from the millimetre flux, using the following equation: 



Mdust 



^vBv(rdust) 



(1) 



where ky = ko(v/vo)^ is the mass absorption coefficient, /3 
parametrises the frequency dependence of ky and Sy is the ob- 
served flux. D is the distance to the source, Tdust is the dust tem- 
perature and Bv(rdust) is the Planck function. 

In order to be consistent with C05, we compute 
the dust masses assuming j3-l, k^-lcm^g^^ at 1.3 mm 
Beckwifli et al. 1990l l and rdust=40 K. 

The uncertainty in the dust mass is computed using error 
propagation. In the following, we discuss the principal system- 
atic errors which are present in the final dust mass uncertainty. 
The mass absorption coefficient is poorly constrained. We 
consider as limit the commonly used value for debris discs 
of ko-1.7cm^g^^ at 800 jum adopted by Zuckerman & 
Becklin (1993), which corresponds to ko=l.Ocm^g^^ at 1.3 mm. 
The mass absorption coefficient alone introduces at least a fac- 
tor 2 in the final uncertainty associated to the mass. In order 
to be consistent with the previous analysis of C05, we adopted 
a temperature of 40 K. This represents a compromise between 
the temperature associated with cold dust in protoplanetary discs 
(r~20-30 K, e.g. Beckwith et al. 1990J and warm dust in debris 
discs (r~40-100 K; Zuckerman & Song |2004al l. 
Since the mass is proportional to Tdusf temperature alone in- 
troduces a factor of 2 in the uncertainty of the masses. 
Another source of uncertainty in the final mass error is intro- 
duced by the source's distance. Nevertheless, the parallax of 
most of the sources older than 20 Myr (which will be in the 
following analysed), has been measured by Hipparcos, and the 
distance uncertainty is directly computed from the parallax un- 
certainty. This contributes an uncertainty of at least ~10% of the 
dust masses. 

The 3cr upper limits to the dust mass have been computed when 
the sources were not detected with a signal to noise ratio of at 
least three. 

The mass of the newly observed sources younger than 20 Myr is 
also plotted as reference Q. Taking into account the uncertainty 
on distance, temperature, measured flux and mass absorption co- 
efficient, the final uncertainty is ~55% of the dust mass. 

3.1. Evolution of the dust mass 

Mass detections and upper limits are presented in Fig. [T] as a 
function of age. In total, we plotted 132 upper limits and 9 detec- 
tions computed from the millimetre measurements. 17 sources 
have been observed in both surveys, C05 and with IRAM (con- 
nected by a dotted line in Fig. [T]i. Due to the higher sensitivity 
of the IRAM bolometer compared to SEST and OVRO, we de- 
tected 4 discs with dust masses down to 5 x 1O"'M0, about one 
order of magnitude smaller than those in C05. 
The other variable which can strongly affect our analysis is the 
age determination of the discs older than 20 Myr. Their ages 
derived using different techniques are summarised in Table |2] 
The final ages are computed as the average and the standard 
deviation is adopted as uncertainty (Table |2|. A more detailed 
discussion of the ages of the whole FEPS sample will be pre- 
sented by Hillenbrand et al. (in preparation, see also Mamajek 



' Note however, that the assumption of optically thin emission at mil- 
limetre wavelength cannot be generalised to all primordial discs 
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Fig. 1. Disc mass versus age of the system. Filled circles are the masses determined from the IRAM detections (this work) and from 
the detections presented by COS. Upper limits are shown as filled triangles (new IRAM observations), empty and filled triangles 
(objects measured by COS and re-observed here with higher sensitivity) and down arrows (from COS). The sources observed in both 
studies are connected by dashed lines. The vertical line separates the sources younger than 20 Myr and evolved systems older than 
20 Myr Primordial discs are shown in this plot only as comparison with systems older than 20 Myr but they are not taken into 
account in any statistical analysis in this paper. 



& Hillenbrand ( l2008l l). 

In the following subsection we present the statistical analysis 
aimed at constraining the physical characteristics of the dust in 
circumstellar discs at ages older than 20 Myr, in order to un- 
derstand the timescale and the dissipation mechanism in debris 
discs. 



Table 2. Ages (in logarithmic scale) of the debris discs detected 
at the millimetre wavelength, with IRAM (this work), SEST and 
OVRO (from COS). They have been derived using different tech- 
niques: the strength of the Call H and K (^^^) and Lithium (Li) 
emission. X-ray activity (X-ray), and the rotational period-age 
(Prot) relation (L.Hillenbrand, private communication). The ages 
are averaged and the standard deviation is adopted as error. 



Metiiod 


HD 377 


HD 8907 


HD 104860 


HD 107146 




8.34 


8.78 


8.44 


8.09 


V sini 


7.99 


8.1 




8.89 


Plot 






7.62 




X-rav 


7.9 


7.88 


8.06 




Li 


7.8 


9.28 


8.08 


8.1 


average 


8.01 


8.51 


8.05 


8.36 


stdev 


0.23 


0.64 


0.34 


0.46 



3.1 .1 . Method: statistical analysis 

A quantitative analysis of the dust mass-age relation was per- 
formed for systems older than ~20 Myr. We carried out a sur- 
vival analysis (Feigelson & Nelson'198S) in order to include the 
information contained in the dust mass upper limits derived from 
the millimetre non-detections. We use the program ASURV Rev. 



1.2 (Lavalley et al. 1992), which implements the survival analy- 
sis methods presented in Feigelson & Nelson ( I198SI ) and Isobe 
et al. OUg. 

The generalised Kendall's tau correlation (Brown et al. 119741 ) 
was used to measure the degree of correlation between age and 
disc mass among debris discs. Such a method includes the anal- 
ysis of both detections and upper limits. The probability that the 
disc mass is correlated with the age is 76%. 
To determine whether there is a different dust mass distribution 
in debris disc systems older than 20 Myr, we used two-sample 
tests. We applied three different tests: the Gehan, logrank, and 
Peto-Prentice tests (e.g. Feigelson & Nelson [T98SI I. In order to 
investigate the evolution of the debris dust we divided the sam- 
ple of objects older than 20 Myr in two sub-samples with al- 
most the same number of objects: the first with ages between 20 
and ~180 Myr, and the second sub-sample with ages between 
-180 MyrandSGyr. 

The three two-sample tests give a probability between 70 and 
83% that younger and older debris systems belong to different 
parent populations in terms of dust mass. The relatively low sig- 
nificance of our results comes from the low number of detected 
sources, compared to the upper limits, and the small range in age 
of the detected debris discs. This leads to a high probability to 
assign a detected disc to a different sub-sample in age. We could 
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not perform any conclusive linear regression to the data points 
since the detected objects span a very narrow range in age and 
masses, and in general there are too many upper limits compared 
to detections. 



3.2. Discussion 

The detection limit of our survey does not allow to draw any firm 
conclusion about the evolution of the debris dust around solar- 
like stars as a function of age. Even more sensitive millimetre 
observations are needed to find statistically reliable trends. Our 
sample is biased toward sources with infrared excess. However, 
with the current sensitivity limit, additional observations of 
sources without infrared excess, which have a smaller proba- 
bility to be detected in the millimetre, would only increase the 
number of upper limits. This would not change the result of our 
statistical analysis. 

We recall that all the mass derivations are valid only if the mass 
absorption coefficient k, which is inversely proportional to the 
mass of the dust in the disc, remains constant with time. 
While in the protoplanetary disc phase we would expect a strong 
change of the mass absorption coefficient with time, in the de- 
bris disc phase we expect an equilibrium phase which is char- 
acterised by an almost constant second generation of dust. The 
assumption of a roughly time-independent dust opacity seems to 
be reasonable as long as we consider discs in the debris phase. 

4. SED analysis of debris discs detected at 350 /im 
and/or 1.2 mm 

The debris discs detected by our surveys offer the opportunity to 
characterise and compare the disc properties around solar type 
stars. This is the case for the following objects: HD 104860, 
HD 8907, HD 377, HD 107146, HD 61005 and HD 191089. 
We compiled the spectral energy distributions from the infrared 
data available in the literature and our new detections at sub- 
millimetre and/or millimetre wavelengths. Synthetic photomet- 
ric points at 13 and 33 fim, derived from the IRS spectrum 
(Hillenbrand et al. 12008 ), were also included in the SEDs. These 
data points allow us to distinguish the infrared excess of the de- 
bris disc from the stellar photosphere, determined by a Kurucz 
model. The IRS spectrum is overplotted in Figures [JE] but was 
not considered during the SED modelling. 
To further characterise the debris discs detected during our sur- 
vey, we compute: 

- the j6 index from their SED to derive the dust grain size present 
in the disc, 

- the collisional and Poynting-Robertson drag timescales 

- the blowout grain size (flbiow) 

- best fit model parameters by modelling the SED. 

jS Index - In the millimetre wavelength range, where 
hv/KT << 1, the Planck function can be approximated with the 
Rayleigh-Jeans relation (By ~ Ikv^Tc^^). The disc is assumed 
to be optically thin with the optical depth Ty{r) - ky 2(r) < 

1. Under these conditions Fy ~ kyV^^ J^'"' l,(r)T(r)rdr . At 
wavelengths A > 0.1mm the mass absorption coefficient ky is 
ky oc where /3 is the dust opacity index. Since the slope of 
the SED in the millimetre wavelength range is proportional to 
V, the measure of the spectral index a enables us to directly 
constrain the /? index via the relation f5 - a - 2 and to derive 
the index The index /? is 2 for interstellar medium grain sizes, 
between and 2 for pebbles of the order of I mm in size, and 
for even larger grains (e.g. Beckwith et al. I2000I I. The slope a 



of the SED in the sub-millimetre/millimetre wavelength range is 
derived, using a linear regression that gives the errors in slope of 
the best-fit. The values for the /3 index of the power law of the 
dust absorption coefficient are always < 1 (see Table [3]l, which 
suggests the presence of dust with at least millimetre size parti- 
cles (e.g. Henning & Stognienko |I9961 Rodmann et al. 120061) . 

Blowout grain size - Small grains are affected by the inter- 
action with the stellar radiation field. This causes a force acting 
on the particles which is parametrized by the ratio of the radia- 
tion force to the stellar gravity. We compute the blowout grain 
size (abiow) in the debris disc systems using the equations pre- 
sented in Hillenbrand et al. (2008, based on Burns et al. 119791 
and Backman & Paresce[1993j): 



flbiow = 0.52jum 



l.lg/cm^ 1 + p 



1.1 (M*/Mo)(r*/5780) 



(2) 



where, p is the grain density, p the grain albedo (0.1 for silicate 
grains) and a the grain radius. Sub-micron sized particles are 
therefore blown out from the disc near a Sun-like star. 

Poynting-Robertson and collisional timescales - Defining 
rji) - fpR/fcoii (Wyatt et al.' 1999 2005 ), where fpR and fcou are the 
Poynting-Robertson drag (fpR) and collisional timescales respec- 
tively, it is possible to distinguish between discs in collisional 
regime (t/q » 1) and Poynting-Robertson drag regime (770 <k 1). 
When rjo ^ \ the disc is dense and the collisions occur faster 
than the P-R drag: the dynamically bound dust remains at the 
same radial location as the planetesimals, while unbound grains 
are blown out and their surface density distribution falls off as 
r"' (where r is the distance from the star). 
The Poynting-Robertson drag (fpR) and collisional timescales 
(fcou) of the debris discs are computed, using the equations pre- 
sented in Hillenbrand et al. (i2008i based on Burns et al. 119791 
and Backman & Paresce[T993): 



fpR = llQyr 



(p/g/cm^)(a/ium)(r/AUY 
(LJLq)(1+p) 



(3) 



where, p is the grain density, p the grain albedo (again 0. 1 for 
silicate grains) and a the grain radius. We assume silicate grains 
of 10 pm in size and a grain density of 2.7 g/cm^ at the inner 
radius of the disc (computed by our model, see Tabled. 
Under the assumption of circular orbits and completely destruc- 
tive collisions between grains of the same size, the collisional 
timescale is computed using the equation: 



fcoU 



1 



f-f- 

\AUJ 9a-{r) sImTIm^ 



yr 



(4) 



where cr{r) is the face-on fractional surface density: for a con- 
stant surface density, cr(r) - 2///n(7?ouT/^in), with R^a the inner 
disc boundary, /?ouT the outer disc radius and / - Ldust/LQ (e.g., 
Backman 2004). 

The collisional timescale has been computed for grains of 10 pm 
in size. According to Wyatt (2005), the ratio fcoii/fpR is propor- 
tional to 1/ Vfl, where a is the grain size. 
These values are summarised in Table |3] 

SED modelling - The SEDs were modelled using the radia- 
tive transfer model of Wolf & Hillenbrand (2003). 
The model assumes that the dust grains are compact spherical 
particles heated only by the direct stellar radiation. The interac- 
tion between stellar radiation and dust particles can be described 
by the following radiative processes: scattering (neglecting mul- 
tiple scattering), absorption and re-emission of stellar radiation 
by dust grains. The disc is then checked to be optically thin to the 
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Table 3. Blowout grain size, P-R drag timescale assuming silicate grains of 10 fim, grain density of 2.5 g/cm^ at a distance of the 
inner radius of the disc (computed by our model, see Table HJ. The yS-index discs only for the stars with the SED sampled in the 
submillimetre/millimetre range (details in section 4) is reported. 





HD 104860 


HD 8907 


HD 377 


HD 107146 


HD 61005 


HD 191089 


abiow.Silj"m] 


0.64 


0.81 


0.54 


0.49 


0.34 


1.14 


fp-R,Si[Myr] 


0.39 


0.5 


0.03 


0.08 


9.85 


1.86 


tcou [Myr] 


0.007 


0.007 


0.003 


0.002 


0.009 


0.005 




0.5±0.7 


-0.1+1.3 




0.4±0.2 







Table 4. Model results, reduced ik^wd assuming a fixed outer radius (^qut) of 150 AU and a maximum grain size (flmax) of 
3 mm. The parameters errors are estimated by the Levenberg-Marquardt algorithm. 



HD 104860 HD 8907 HD 377 HD 107146 HD 61005 HD 191089 

Md„„/Me 0.082+0.007 0.040+0.005 0.058+0.013 0.110+0.008 0.261+0.023 0.110+0.037 

R,n [AU] 21+5 27.8+13.9 6.1+0.8 10.2+0.8 95.6+23 79+98 

flniin [//m] 8+2 6+2 14+5 8.6+1.2 0.4+0.3 0.23+0.78 

2.6 4.3 5.0 4.5 16.3 8.7 



stellar radiation and to the dust re-emission at all wavelengths. 
The contribution to the SED from a single dust grain results from 
the integration of the scattering and re-emission processes over 
all wavelengths. Finally the emergent spectrum is the sum of all 
the dust grain contributions. 

The model parameters are the radiation emitted by the central 
star, the total disc mass, the disc size, the radial density dis- 
tribution, the minimum and maximum grain size in the disc, 
the grain size distribution, and the chemical composition of the 
grains. The radial density distribution is described by a power- 
law, n(r) oc r^i where q - I corresponds to a disc with a constant 
surface density E(r) oc r° . The grain size distribution is also de- 
scribed by a power-law n{a) oc a^P where the canonical value 
p = 3.5 characterises a size distribution initially produced by a 
coUisional cascade. We always use this grain size distribution in 
our models. For each grain size a temperature distribution over 
the disc size is computed. 

The model parameter space has been extensively analysed by 
Wolf & Hillenbrand ( I2003I I. They found that the increase of the 
inner radius Ri„ causes a loss of the warm dust that is mainly re- 
sponsible for the near-infrared (NIR)/ mid-infrared (MIR) shape 
of the SED; the excess emission over the stellar photosphere 
starts at longer wavelengths. This is more pronounced when the 
minimum grain size a,ni„ increases, because the main contribu- 
tion to the NIR/MIR spectrum at high temperatures comes from 
small grains. The presence of an inner gap causes the loss of 
warm dust which is mainly responsible for the NIR/MIR shape 
of the SED. Increasing the size of the gap, the MIR flux de- 
creases and the excess is shifted to longer wavelengths. Keeping 
the disc mass constant, the flux in the millimetre region increases 
slightly, but the net flux is smaller compared to a disc without an 
inner gap. This is because the fraction of the stellar flux absorbed 
by a single dust grain decreases with increasing radial distance 
from the star. 

We assumed astronomical silicate grains (optical data from 
Draine & Lee |1984| and Weingartner & Draine 2001) since they 
are expected to be the main component of the dust in protoplan- 
etary discs (Pollack et al. |1994| ). 

Since the outer radius cannot be constrained from the SED alone, 
it was assumed to be 150 AU, consistent with the location of the 
most distant object of the Kuiper belt in our Solar System (e.g. 
Gladman et al. 120021 ). The choice of such a value is supported by 



a recent resolved millimetre map of one of the debris discs anal- 
ysed (Corder et al. 120091 ). where the millimetre emission is ex- 
tended up to 150 AU from the central star The maximum grains 
size is assumed to be 3 mm, since the /3 index measurements 
suggested the presence of at least millimetre-sized particles in 
the discs; this is the grain size which strongly contributes to the 
millimetre emission that we detect. The Levenberg-Marquardt 
method was applied to find the best-fit parameters of the model 
leaving the inner radius (Rm), minimum grain size fflminj ond 
dust mass (Mdustj as free parameters. 

The reduced associated with each best-fit model, was com- 
puted for each fit normalising the standard statistic to the 
number of free parameters of the fit (Bevington |1969l l. These 
parameter uncertainties are estimated from the covariance ma- 
trix (inverse of the curvature function) by the Levenberg- 
Marquardt algorithm. The modelling results are listed in Table|4] 
The degeneracy of the SED models was checked using the cr- 
confidence regions delineated by lines of constant Ax^- The 
model parameter space (^inxMdustXflmin) was explored around 
the best fit values. The 2-dimensional confidence regions in 
Figs. |2ll2] are the projections of the 3-dimensional isosurfaces 
defined by fixed Ax^ (e.g. Press et al. 1992). The shape of the 
confidence regions confirm an intrinsic degeneracy of our model 
probably between inner radius, minimum grain size and dust 
mass. 



4.1. HD 104860 

HD 104860 is an F8 star at a distance of 47 pc from the 
Sun. It was classified as a y Doradus variable candidate by 
Handler (il999i) . Wichmann et al. (120031) identified the object 
as a young zero age main sequence star using Lil in absorption 
(Wichmann et al. 120031) . In Table |2] we summarise the ages de- 
rived from different methods (L. Hillenbrand, private communi- 
cation). The average age is ~140 Myr, the value that we adopted 
in our analysis. 

The SED of HD 104860 is presented in Fig.|2] Infrared-excess 
emission has been detected by Spitzer, suggesting the presence 
of debris material around this main sequence star (Carpenter et 
al. l2008l ). 

Early modelling of the spectral energy distribution of 
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Fig. 2. SED of HD 104860 (upper plot (a)): Detections are 
shown as empty diamonds while 3cr upper limits are shown with 
empty triangles. It includes the Spitzer/IRAC and Spitzer/MIPS 
photometry from 3.6 to 160 //m, synthetic photometry at 13 
and 33 fim obtained from the IRS low resolution spectrum 
(Hillenbrand et al. (I2008I I. our CSO and IRAM detections at 
350 yum and 1.2 m m, SC UBA detections at 450 and 850 yum 
(Najita & Williams "2^05^, and an OVRO 3.1 mm upper limit 
(Carpenter et al. i2005i . The solid line shows the best fit model 
(see details in section 4). The lower plots (b, c, d) represent the 
2D projections of the 3D surfaces of constant Ax^ obtained vary- 
ing at the same time the dust mass, minimum grain size and inner 
disc radius. 



HD 104860 was presented by Najita & Williams ( |2005] ). 
Treating dust grains as simple blackbody radiators, their sub- 
millimetre observations at 450 and 850 fim indicate a dust tem- 
perature of 7'ciust=33 K and a mass of the disc of 0.16 M®. Our 
new detection allows us to sample the SED at longer wave- 
lengths and hence constrain the slope of the far-infrared to mil- 
limetre excess. The best-fit model parameters are reported in 
Table|4]and the computed best-fit model is plotted in Fig.|2] The 
confidence regions of the fit (see lower plots in Fig. |2]i show 
some degree of degeneracy of the model. The dust mass is con- 
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0.02 f ................... 1 
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Mininnum grain radius [/xm] 

Fig. 3. As Fig. |2] but for HD 8907. The I SO photometry at 
60 and 90 jjm are from Silverstone (l2000l l. the SCUBA de- 
tections at 450 and 850 /im from Najita & Williams (2005|, 
IRAM detection at 1.2 mm from this work and the OVRO 
2.6 and 3.1 mm upper limits from C05. The IRAS detections 
JHelou & Walker 19881 ). shown by empty squares, are systemat- 
ically larger than the Spitzer photometry, due to the bigger beam 
and they were not taken into account during the model fitting. 



strained between 0.07 and 0.15 M®, while the inner radius 
is located between 10 and 50 AU and the minimum grain size 
is between 5 and 15 fxm. The dust mass computed by the model 
(see Table |4| is of the same order of magnitude as the masses 
computed using a simple black-body model (Najita & Williams 
^2005 ) and Eq. [T] The minimum grain size obtained is one order 
of magnitude higher than the blowout grain size of the system. 
The confidence regions delineated by constant Ax^ contours sug- 
gest that a minimum grain size of the same order of magnitude 
of the blowout size is a poor fit to the data. 

4.2. HD 8907 

HD 8907 is a F8 star, located at 34 pc from the Sun, which has 
already reached the main sequence. In our analysis we adopt an 
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Fig. 4. As Fig. [2] but for HD 377. In the sub- 
millimetre/millimetre wavelenght range HD 377 was only 
detected by IRAM at 1.2 mm (this work). The OVRO 2.6 and 
3.1 mm upper limits are from COS. 



age of the system of 680 Myr, which is the average of the ages 
derived from different methods (Table |2l L. Hillenbrand, private 
communication). 

The SED of HD 8907 is presented in Fig. [3] The infrared excess 
detected by ISO (Zuckerman & Song |2004al l suggested the pres- 
ence of debris dust. Under the assumption that the dust particles 
radiate as a black body, a temperature of ~ 60 K was inferred for 
particles orbiting at ~30 AU from the central star (Zuckerman 
& Song .2004 a). The SED infrared excess has also been detected 
by Spitzer between 3.6 and 70 fim by Kim et al. ( 120051 ), and they 
found an inner radius of ~48 AU with a temperature of ~ 48 K. 
Our new detection allows us to sample the SED at longer wave- 
lengths and hence constrain the slope of the far-infrared to mil- 
limetre excess. The best-fit model parameters are reported in 
Tableland the computed best-fit model is plotted in Fig. [3] The 
confidence regions of the fit which show some degree of degen- 
eracy of the model (see lower plots in Fig. [3j: the dust mass is 
constrained between 0.03 M® and 0. 1 M®, while the inner radius 
is located between 15 and 30 AU and the minimum grain size is 
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Fig. 5. As Fig. |2]but for HD 107146. The source was not ob- 
served during our survey with IRAM, but it was recently de- 
tected by CARMA at 1.3 m m (Corder et al. [20091 . The IRAS 
detections dHelou & Walker 1988) , shown by empty squares, are 
systematically higher than the Spitzer photometry, due to the 
bigger beam and they were not taken into account during the 
model fitting. 



between 5 and 10 /urn. 

As in the case of HD 104860, the dust mass computed by the 
model of HD 8907 (see Table|4]i is comparable to the value com- 
puted using Eq. [T] The minimum grain size is between 4 and 
8 yum. This is one order of magnitude higher than the blowout 
grain size of the system and is consistent with previous mod- 
elling (Kim et al. 12005b . The inner radius computed by our 
model is poorly constrained since it has an uncertainty of 50%. 

4.3. HD 377 

HD 377 is a G2V star at a distance of 40 pc which has al- 
ready reached the main sequence. It is surrounded by a de- 
bris disc with a high fractional luminosity (LiR/Lstar), 4.0x10"'* 
(IMoor et al. 20061 1. The ages of HD 377 derived from diff'er- 
ent methods (L. Hillenbrand, private communication) are sum- 
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Fig. 6. As Fig. El but for HD 61005. In the sub- 
millimetre/millimetre wavelenght range HD 61005 was only de- 
tected by CSO at 350 jum. SEST upper limit at 1.2 mm is from 
COS. 



marised in Table |2] The average age is ~1 10 Myr, the value that 
we adopted in our analysis. 

Pascucci et al. ( I2006I I investigated the presence of gas in the disc 
of HD 377, but they did not detect any gas emission lines in 
the infrared nor in the millimetre wavelength range. The SED 
was compiled from the literature including our new detection at 
1.2 mm (Fig. |4]i. The best-fit model parameters are reported in 
Tableland the computed best-fit model is plotted in Fig.|4] The 
confidence regions of the fit (see lower plots in Fig. |4]i show a 
high degree of degeneracy of the model compared to the other 
debris discs modelled. This results from the poor constraints on 
the SED in the sub-millimetre region. 

As for HD 104860 and HD 8907, the dust mass computed by the 
model (see Table |4]l is comparable to the value computed using 
Eq.[T] The inner radius is the smallest in our sample, but the disc 
seems to lack small dust grains with sizes up to 14 fim. The min- 
imum grain size obtained is one order of magnitude higher than 
the blowout grain size of the system. 
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Fig. 7. As Fig. |2]but for HD191089. CSO detection at 350 //m 
and IRAM upper limit at 1.2 mm are from this work. OVRO 
3.1 mm upper limit is from Carpenter et al. (l2005i l. The IRAS 
detections (IHelou & Walker 198 8 J , shown by empty squares, are 
systematically higher than the Spitzer photometry, due to the 
bigger beam and they were not taken into account during the 
model fitting. Since there is no IRS spectrum of the source con- 
straining the beginning of the infrared excess close to the stellar 
photosphere, the IRAS detection at 12 yum was taken into ac- 
count during the SED fitting. 



4.4. HD 107146 

HD 107146 is a G2 V star located 28.5 pc from the Sun. 
The star is a candidate periodic V-band photometric variable 
( |Koen & Eyer 20021 ). From the average between the ages derived 
with different methods (Table|2]l, we obtain an age of ~230 Myr 
HD 107146 is also the first debris disc around a solar-type 
star resolved in scattered light by the Hubble Space Telescope 
(lArdila et al. 2004l i. They resolved a ring-like disc with most 
of the material concentrated between 86 and 130 AU, using 
a coronograph of 1.8"in radius. The disc was also marginally 
resolved at 450 and 850 fim and the observations suggested 
the presence of an inner hole with a radius larger than 3 1 AU 
( IWilliams et al. 20041 ). Such an inner hole would also explain the 
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lack of an IRAS 25 fim excess and a far-infrared excess at 60 and 
100 iim (Silverstone |2000] l. COS presented OVRO 3 mm images 
of this source . HD 1 07 1 46 has also been resolved during our sur- 
vey at 350 yum with CSO and at 1 .3 mm with CARMA (Corder et 
al. 120091 1. The millimetre image shows a clumpy disc extended 
between ~60 AU and ~150 AU with a symmetric peak of the 
emission at ~91 AU from the central star The SED was com- 
piled from the literature, including our new detection at 350 fim 
(Fig. Is). HD 107146 is the only disc detected at 3.1 mm. The 
best-fit model parameters are reported in Table |4] and the com- 
puted best-fit model is plotted in Fig. |5] The degree of degen- 
eracy of the model is shown in the confidence regions of the fit 
(see lower plots in Fig.|5]l. The dust mass is constrained between 
0.08 Me and 0.15 M®, while the inner radius is located between 7 
and 15 AU and the minimum grain size is between 8 and 15 fim. 
It is interesting to notice that the inner radius computed by our 
model, is more than 7 times smaller than the inner radius of the 
ring-like disc in the HST images. This is because different grain 
populations are traced by our model and the HST image: while 
the HST traces only the distribution of the sub-micron size pop- 
ulation (below flbiow), our model includes the re-emission from 
grains from micron to millimetre in size. This explanation is also 
supported by the resolved images in the sub-millimeter wave- 
length range of HD 107146 which suggest the presence of grains 
at radii smaller compared to the ring-like structure seen in the 
HST image. 

4.5. HD 61005 

HD 61005 is a G3/G5V star with a Hipparcos distance of 
35+1 pc, located in the local bubble, a region that is thought 
to be almost completely free of diffuse dust clouds (Franco et al. 
|1990 i From an average of the chromospheric activity available 
from the literature, and a new calibration of the chromospheric 
activity-age relation (FEPS, private communication), we adopt 
an age of 135 Myr. The SED was previously modelled by Hines 
et al. ( 120071 ): under the assumption that dust particles of ~10 fim 
radiate as a black body with a temperature of ~50-70 K, they 
obtained a minimum distance for the circumstellar material of 
~7 AU. This star was also monitored during a radial velocity 
survey looking for planets around active stars, but no planets 
were detected (Setiawan et al. 2008). HST/NICMOS observa- 
tions (Hines et al. 2007 ) reveal dust-scattered starlight extend- 
ing to distances of ~240 AU from the occulted star The struc- 
ture is strongly asymmetric about its major axis, but is mirror- 
symmetric about its minor axis; morphologically, the object re- 
sembles a wing-spread moth with the star as the head. 
HST scattered light images reveal a swept shape of the disc with 
an inner radius < 1 AU and an outer radius up to 240 AU (Hines 
et al. 120071 ). Such emission is thought to be associated to the lo- 
cal interstellar medium that scatters the stellar light, so that the 
movement of HD 61005 through the medium causes the swept 
shape of the disc (Hines et al. |20071 l. These observations give a 
complementary picture of the circumstellar material present in 
the disc since different grain populations are traced in our model 
and in the scattered light images. The SED was compiled from 
the literature including our new detection at 350 fim (Fig. |6]). 
The best-fit model parameters are reported in Table |4] and the 
computed best-fit model is plotted in Fig.|6] The confidence re- 
gions of the fit (see lower plots in Fig. |6]l show that the best 
fit parameters represent only a local minimum of the model pa- 
rameters space. This is probably due to the fact that the model 
assumptions are not valid in this case. As suggested from the 
HST images, the object seems to be located in the local interstel- 



lar medium material, which probably introduces an additional 
small grain population which contributes to the scattered light in 
the SED. Our model assumes that all the detected flux is stellar 
light scattered and/or absorbed and re-emitted only by the dust 
particles in the disc. 

4.6. HD 191089 

HD 191089 is an F5V star with an Hipparcos distance of 
54+3 pc from the Sun. This star was selected for the FEPS pro- 
gram because it was known to have a debris disc based on ISO 
and IRAS measurements and it was not considered in our sta- 
tistical analysis presented in section 3. An upper limit of the 
age (~3 Gyr) for the system was presented by Nordstrom et 
al. (2004), using main-sequence isochrones. We adopt an age 
of -300 Myr (Hillenbrand et al. 2008 ). The SED compiled from 
the literature and the Spitzer spectro-photometry (Carpenter et 
al. 2008) is presented in Fig. [T] The best-fit model parameters 
are reported in Table |4] and the computed best-fit model is plot- 
ted in Fig.|2] The confidence regions of the fit show some degree 
of degeneracy of the model (see lower plots in Fig.|2l). We note 
that the model result is not significant at all. This is due to the 
few data points which constrain the SED. 

4. 7. Modelling discussion 

Comparing the collisional timescale to the Poynting-Robertson 
timescale, we conclude that all the detected debris disc are 
systems in the collision dominated regime. 
The SED of the debris discs has been modelled leaving 3 
parameters simultaneously free: inner radius, dust mass and 
minimum grain size. 

We notice that the slope of the SEDs seems to change from the 
sub-millimetre to the millimetre wavelength range, but is well 
reproduced by our models within the errors. Including grains 
larger than 3 mm in the model would decrease the steepness of 
the slope in the sub-millimetre, but would not reproduce at the 
same time the sub-millimetre and millimetre fluxes. 
A degeneracy between the inner radius and minimum grain 
size can still exist: in this case the cr-contours plots would 
represent only a local minimum in the parameter space, instead 
of the absolute solution. Only a resolved image of the entire 
grain population which delineates the inner radius of the disc 
will allow to break the degeneracy between smallest grain size 
present in the disc and inner radius. 

We now discuss how well constrained the fixed model parame- 
ters are (i.e. Rom, flmax, dust composition and size distribution) 
and in particular if and how they can alter our conclusions 
about Mdust, ^in ™d We are not going to further discuss 
the cases of HD 61005 and HD 191089 because our modelling 
was not significant at all, nor the case of HD 107146 because 
the assumption on the outer radius is supported by millimetre 
resolved images (Corder et al. 120091 ). 

- The outer radius was assumed to be 150 AU. While in the 
case of HD107146 such an assumption is supported by resolved 
millimetre images, in the other cases, we need to investigate the 
robustness of our model results. While the most distant objects 
of the Kuiper belt have been observed up to 150 AU from the 
Sun (e.g. Gladman et al. 120021 ). the edge of the Kuiper belt is 
commonly considered close to the 2:1 resonance with Neptune 
at about 50 AU from the Sun (e.g. TrujiUo & Brown 20011. 
There is also no reason why all discs would be expected to have 
the same outer edge. For these reasons we use five different 
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Table 5. Model results leaving as free parameters using Mdust/Mg, Ri„ and Omin for 6 different outer radii: 1000, 500, 300, 150, 100 
and 50 AU. 



HD 104860 



Rout 


1000 AU 


500 AU 


300 AU 


150 AU 


100 AU 


50 AU 




5.843±0.652 


2.930±0.245 


1.774 0.133 


0.082±0.007 


0.053±0.005 


0.039±0.004 


^in [AU] 


185± 23 


20.3±7.8 


20±7 


21±5 


22±2 


30±3 




0.07±0.03 


4.5±0.9 


6±2 


8±5 


11±2 


21±4 




6.1 


7.0 


4.1 


2.6 


2.7 


3.0 


HD 8907 


Rout 


1000 AU 


500 AU 


300 AU 


150 AU 


100 AU 


60 AU 




3.097±184 


1.614±0.099 


0.931±0.050 


0.040±0.005 


0.024±0.002 


0.017±0.003 


R,n [AU] 


89±14 


29±14 


19.4±5.4 


27±17 


27±7 


41±15 


amin [/Jm] 


0.5±0.08 


3.1±0.8 


5.6+1.1 


6±2 


8±2 


7.7±2.6 




22.8 


17.8 


9.2 


4.3 


3.6 


3.5 


HD377 


Rout 


1000 AU 


500 AU 


300 AU 


150 AU 


100 AU 


60 AU 


Md„s,/Me 


3.655±0.507 


2.649±0.322 


1.339 ±0.367 


0.058±0.013 


0.034±0.006 


0.039±0.012 


^in [AU] 


70±12 


4.6±1.1 


4.6±0.65 


6.1±0.8 


38±8 


51±8 


Onun [/Jm] 


0.02±0.001 


6.5±1.5 


9.1±3.5 


14±5 


14±4 


38±20 


•^red. 


8.6 


7.7 


4.6 


5.0 


395. 


395. 



outer radii (1000, 500, 300, 100 and 50 AU, or 60 AU when 
the fitting routine was not able to find any reasonable inner 
radius smaller than 50 AU), leaving free the dust mass (Mjust), 
the inner radius (Rin) and the minimum grain size a^in. The 
results are summarised in Table |5] We notice that, in the case 
of HD 104860 and HD 8907, changing the outer radius from 
500 to 50 AU affects only the dust mass. In the case of HD 377 
fixing the outer radius to 100 and 50 AU the model fitting is no 
longer satisfactory, while we find a better result increasing the 
outer radius up to 500 AU (see different in Table |5]). The 
minimum grain size remains at least one order of magnitude 
bigger than the blowout size, despite the case of an outer radius 
of 1000 AU where it becomes of the same order of magnitude 
or smaller than the blowout size. 

- The maximum grain size has been fixed to 3 mm; despite the 
fact that a collisional cascade generates a distribution of grain 
size which does not stop at the millimetre size, we include in our 
model the maximum particle size which contributes significantly 
to the emission which we detect during our observations. We 
remind the reader that the mass derived by our model represents 
the mass of the grains smaller than 3 mm present in the disc and 
that increasing the maximum grain size would increase the disc 
mass. 

- The dust has been assumed to be made up of astronomical 
silicate; this is expected to be the main component in pro- 
toplanetary discs (Pollack et al. |1994| l. Although a different 
composition of the dust can be present in the disc, Spitzer 
infrared spectroscopic observations do not allow to constrain 
the material properties due to the lack of diagnostic features, 
caused by large grain sizes (see the IRS spectra overplotted in 
the SEDs, Fig. 2-7). A different dust composition would change 
the blow-out size as well as the minimum grain size present in 
the disk computed by our model. A detailed investigation of 
different dust compositions is beyond the goal of this paper. 

- We finally recall that a more realistic grain size distribution 
can also affect our modelling results. Departures from the 
canonical size distribution we assumed, have been highlighted 



by collisional evolution studies of spatially resolved debris disks 
(e.g. Thebault & Augereau |2007l l. 

Despite previously discussed degeneracies, the objects 
HD 104860, HD 8907 and HD 107146 show the inner 
part of their discs cleared from small micron-sized dust grains 
because of the lack of near/mid-infrared emission (Wolf & 
Hillenbrand I2003I I . These inner holes can be maintained by the 
presence of a planet-sized body (Roques et al. 1994), which 
avoids the filling of the central cavity by Poynting-Robertson 
drag. 

5. Summary and conclusions 

We have carried out two deep surveys at 350 yum and 1.2 mm 
of circumstellar discs around solar-type stars at ages between 
3 Myr and 3 Gyr. These new observations are an order of magni- 
tude more sensitive to dust emission than previous observations 
of C05. The dust disc masses have been computed from the mil- 
limetre emission, where the discs are assumed to be optically 
thin. A survival analysis of the dust disc masses as a function of 
time has been carried out of systems older than 20 Myr, includ- 
ing the dust mass upper limits. The spectral energy distributions 
of the debris discs detected by our sub-millimetre and millime- 
tre surveys have been modelled. We draw the following main 
conclusions from our work: 

1. The Kendall's tau correlation yields a probability of 76% 
that the mass of debris discs and their age are correlated. 
Similarly, the three different two-sample tests gives a proba- 
bility between 70 and 83% that younger and older debris sys- 
tems belong to different parent populations in terms of dust 
mass. Our result on the relation between dust mass and age 
is limited by the sensitivity of our millimetre survey. Deeper 
millimetre observations are needed to confirm the evolution 
of debris material around solar-like stars. 

2. The spectral energy distributions of the debris discs detected 
at 350 /vm and/or 1.2 mm were modelled. We found a de- 
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generacy in the best fit parameters which will only be bro- 
ken with high-resolution images which resolve the entire 
disc. Nevertheless, this approach allows us to identify de- 
bris discs with an inner region that has been evacuated from 
small micron-sized dust grains. 
3. In the case of the detected debris discs, the comparison be- 
tween collision and Poynting-Robertson timescales suggests 
that the debris discs are collision dominated. 
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Table 1. The stellar properties, e.g. spectral type, temperature, luminosity and age range reported in the first columns are from Meyer et al. (2006)1, 
Hillenbrand et al. ( I2008I I and Meyer et al. (2008 1. CSO 350 yum and IRAM 1.2 mm measured fluxes and upper limits are reported in the sixth and 
seventh column. The detected sources are highlighted in bold. The errors and the upper limits refer to the statistical errors (flux rms) and do not 
include the calibration uncertainty. In the ninth column are the references to the stellar distances. 
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Table 1. continued. 



Source 


Spectral 


Teff 


logiL-klLo) 


log{Age) 


350/im 


s ** 

1.2 mm 


Distance 


Ref. 




Type 


[°K] 






[mJy] 


[mJy] 


[pc] 




HD 38207 
HD 117176 
HD 218738 


F2V 


6769 


0.72 


8.0-8.5 


<21.9 


<0.33 
<14.1 
<3.3 


127+25 


(3) 



*: The calibration error is between 12 and 22% of the flux. **: The calibration error is between 11 and 16% of the flux. Distance references: ' 
de Zeeuw et al. ( 1999); ^Ferryman et al. J 19971 ; ' Hillenbrand et al ( 12008 1 ; Pinsonneault et al. ( |1998t ; ^ based on recent Pleiades distances: de 
Zeeuw et al. a999> , Zwahlen et al. ( 12004) . Soderblom et al. ( |2004| | and Pan et al. t2004< . 



